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Abstract. Inspired by the 'computable cross norm' or 'realignment' criterion, we propose a 
new point of view about the characterization of the states of bipartite quantum systems. We 
consider a Schmidt decomposition of a bipartite density operator. The corresponding Schmidt 
coefficients, or the associated symmetric polynomials, are regarded as quantities that can be 
used to characterize bipartite quantum states. In particular, starting from the realignment 
criterion, a family of necessary conditions for the separability of bipartite quantum states is 
derived. We conjecture that these conditions, which are weaker than the parent criterion, can 
be strengthened in such a way to obtain a new family of criteria that are independent of the 
original one. This conjecture is supported by numerical examples for the low dimensional 
cases. These ideas can be applied to the study of quantum channels, leading to a relation 
between the rate of contraction of a map and its ability to preserve entanglement. 



PACS numbers: 03.67.-a, 03.67.Mn 



1. Introduction 



The relation between the state of a composite quantum system as a whole and the 
configuration of its parts is a very peculiar feature of quantum theory. As recognized since 
the early stages of development of the theory [Q3 121 El, this is a consequence of the tensor 
product structure of the state space of a composite quantum system. This feature of quantum 
mechanics has its most evident manifestation in the fact that it allows the presence of non- 
classical correlations, i.e. of entanglement, between the subsystems of a composite system. 
Nowadays, we may say that quantum entanglement is not only regarded as a key for the 
interpretation of quantum mechanics or as a mere scientific curiosity, but also as a fundamental 
resource for quantum information, communication and computation tasks [@J Q. However, 
despite the great efforts made by the scientific community in the past decades, there are still 
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several open issues regarding the mathematical characterization of composite quantum states, 
even in the 'elementary case' of a bipartite system with a finite number of levels. 

A major challenge is to characterize those states of a bipartite system that are entangled. 
According to the definition due to R. F. Werner J6]|, entangled (mixed) states differ from 
separable states since they cannot be prepared, not even in principle, from product states 
by means of local operations and classical communication only. In mathematical terms, a 
(mixed) state p — a positive (trace class) operator of unit trace — in a composite Hilbert 
space H = H A ® H B is called separable if it can be represented as a convex sum of product 
states: 



with pi > and J2iP* = ^ otherwise, p is said to be entangled. We remark that, if p is 
separable, decomposition (OQ) is in general not unique, and the smallest number of terms in the 
sum (usually called cardinality), due to Caratheodory's theorem, is not larger than the squared 
dimension of the total Hilbert space of the system H (see 0). 

Since quantum entanglement is a very important subject, also in view of its several 
potential applications, separability criteria are regarded as extremely precious tools. Among a 
plethora of proposed separability criteria — i.e. suitable conditions satisfied by all separable 
states whose violation allows to detect entanglement (see, for instance, [[HI [9l [TOl [TTJ ) — the 
present contribution is mainly inspired by the criterion that was proposed in lfT2l with the 
name of 'realignment criterion' (RC) and in lfT3l with the name of 'computable cross norm' 
criterion. As we will try to argue, the RC brings attention to the role played by the Schmidt 
coefficients [HO of a bipartite quantum state in the characterization of entanglement. Trying 
to shed light on this role will be the main goal of our contribution. 

The paper develops along the following lines. In section [2l we introduce a 'Schmidt 
equivalence relation' in the set of states of a bipartite quantum system, and we show the 
link between this notion and the RC. Section [3] is devoted to the characterization of the 
Schmidt equivalence classes. We follow two different approaches: the characterization of 
some groups acting on the Schmidt equivalence classes and the analysis of the local geometry 
of these equivalence classes regarded as manifolds. A family of separability criteria is 
presented in section HI which are extensions of the RC. These criteria are based on the 
'symmetric polynomials' in the Schmidt coefficients, and are weaker than the parent criterion. 
In section [5l the well known correspondence between quantum states and quantum maps 
(see, e.g., |fl4l[T5l[T6ll ). i.e. completely positive trace-preserving (CPT) maps, is considered, 
and a straightforward application of the derived family of separability criteria to the study 
of CPT maps is discussed. Through this correspondence, separable states are associated to 
entanglement breaking (EB) channels ifTTll . The proposed family of criteria, applied to this 
context, leads to a purely geometrical characterization of EB maps. In section[6l we formulate 
the conjecture that the proposed criteria can be strengthened in order to obtain new necessary 
conditions for separability which are independent of the parent RC. Numerical examples in 
support of this thesis are provided in section [7] for low dimensional bipartite systems. 




(1) 
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2. Schmidt equivalence classes of states of a bipartite quantum system 

Let us consider a bipartite, finite-dimensional, complex Hilbert space H = H A <8> H B — 
7Y A ^ <D Na , H b = <C Nb , N a ,N b > 2 — and the corresponding real vector spaces of 
Hermitian operators £ R (7i), £ R (7i A ), and C^(H B ) (the spaces of observables) in H, H A 
and H B , respectively, that are naturally endowed with a scalar product, namely, the bilinear 
Hilbert-Schmidt (HS) product: 

(A, B) us = tr(AB), (A, B Hermitian). (2) 

In particular, the density operators — i.e. the positive operators of unit trace — in H (7i A 
and H B , respectively) can be regarded as elements of the real vector space Ar(W) (£r(H a ) 
and C]n(H B ), respectively), in which they form a convex body that will be denoted by V(H). 
Observe that C U (H) = £ U (H A )®C U (H B ), with £ U (H A ) = H Na ' 2 and C U (H B ) = R N * 2 . The 
HS product allows to write a (nonunique) Schmidt decomposition |[T8l of a density operator 
p G Z>(ft),i.e0 

d 

p = ^A fc F fe A ®F fc B , A!>A 2 >...> A d >0, (3) 
fc=i 

where 

tr(F*F£) = 5 hk = tr(F°F k B ), h, k 6 {1, . . . d}, (4) 

and the real positive numbers {Afc}^ =1 are the (uniquely determined) Schmidt coefficients (in 
short, SC's). Note that the number of terms in the sum equals d = min{A^ A 2 , iV B 2 } (we will 
also set D = max{iV A 2 , iV B 2 }). The definition of the SC's of a bipartite density operator is the 
natural generalization of the standard definition for pure states (see, for instance, |fT9l ). It is 
worth stressing that, since the operators forming the orthonormal systems {F^}'1 =1 , {F fe B }£ =1 
(in £r(7Y a ) and £ R (7i B ), respectively) are Hermitian, the operators {F fc A Cg) F fc B }^ =1 are 
observables. Hence — at least in principle — the SC's are physically measurable quantities: 

\ k = tv(pF£®F k B ). (5) 

The set of 'local' operators { F£ ® F k B } k=1 are also referred to as local orthogonal observables 
EOl . Decomposition © has been recently considered — see [12711 — in connection with the 
formulation of new separability criteria. 

We observe that the convex body V(H) can also be regarded as immersed in the complex 
vector space £(H) of linear operators in H, vector space that can be endowed with the 
sesquilinear HS product (denoted, again, as (•, -)hs)- Then, one can consider a Schmidt 
decomposition of a density operator p in H with respect to the complex Hilbert space C(H). 
It is clear that such a decomposition will contain the same SC's as decomposition ©, but this 
time will involve an orthonormal system of, in general, non-Hermitian operators. 

Given a bipartite density operator, one can uniquely determine its SC's. On the other 
hand, it is clear that the SC's do not identify a unique quantum state. It is then natural to 
formulate the following definition: 

| We remark that actually any operator in (H) admits a Schmidt decomposition. 
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Definition 1 (Schmidt equivalence relation) We say that two bipartite density operators are 
Schmidt equivalent if they share the same set of Schmidt coefficients. 

Then the convex set of density operator in the Hilbert space H — which will be denoted by 
V(H) — is partitioned into Schmidt equivalence classes. 

It is known that a bipartite pure state is completely characterized, with respect to 
entanglement, by the corresponding SC's ||22|| . Although the same characterization cannot be 
extended to a generic state, it is reasonable to suppose that there exists some relation between 
the SC's of a bipartite density operator and the entanglement properties of this state, and 
to address the following questions: "What relevant properties are encoded by the SC's of a 
bipartite density operator, and how one can characterize the Schmidt equivalence classes ? " In 
the following, we will try to analyze these questions and to provide some reasonable answers. 

A first observation is that the SC's determine the purity of a state. Let us recall that the 
purity is defined as the trace of the square of the density operator: 

V{p) := tr (p 2 ) ; (6) 

hence, V(p) G]0, 1], Vp G V(H). It follows from the definition of Schmidt decomposition 
that the purity equals the sum of the squares of the SC's: 

1 d 

jjjj- < P{p) = (p,P>hs = $> 2 < I- (7) 

k=l 

Thus, the purity is the simplest property which is (completely) described by the SC's of a 
density operator. Another relevant fact is the existence of a link between the separability of a 
density operator and its SC's. 

Indeed, the 'realignment criterion' (in short RC see |[T2l [T3l : see also lfT9l . where a 
generalization of the RC is obtained) establishes a necessary condition for the separability of 
a quantum state (or a sufficient condition for the nonseparability). It can be formulated in 
various equivalent ways. From our point of view, it can be regarded as a condition on the SC's 
of a separable density operator p. Precisely, it imposes an upper bound for the sum of its SC's: 



Theorem 1 (the 'realignment criterion') If a bipartite density operator p is separable, then 
its Schmidt coefficients {Xk} A k= i satisfy the following inequality: 

A 

J> fc <l. (8) 

k=l 

The RC is easily implementable. In particular, in [fT2ll it has been introduced the notion of a 
realigned matrix p R associated with the bipartite density operator p in order to compute the 
l.h.s. of inequality ([8]). Fixed orfhonormal bases {|^)} n =i,...Af A and {\y)}v=i,...N B in the local 
Hilbert spaces H A , H B , respectively, and assuming that pi mf j){nv) is the representative matrix 
of the density operator p with respect to the product basis {\n) <g> \v)Yn^i n — wnere ijnp), 
(nu) are double indexes — the corresponding realigned matrix (with respect to the given 
basis) is defined as 

(9) 
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It turns out that the SC's of p are the singular values of the realigned matrix p R , see ||T9ll . 
Precisely, consider a singular value decomposition of the realigned matrix, i.e. 

p R = WAV, (10) 

where hi and V are unitary matrices, belonging respectively to the unitary groups U(iV A 2 ) 
and U(iV B 2 ), and A is a rectangular matrix such that its nonvanishing entries are positive and 
placed along the principal diagonal only. Then, the diagonal entries of A are the singular 
values of p R , hence, the SC's of p. At this point, observing that the singular values of the 
realigned matrix p R coincide with the eigenvalues of the positive matrix |p R | = y / p R V R , one 
concludes that inequality © can also be written as 

||p R || tr = tr(|p R |) < 1. (11) 

Given a density operator p 6 T>(H), we will denote by & p the Schmidt equivalence class 
containing p. Beside the Schmidt equivalence class & p , we will consider the extended Schmidt 
equivalence class containing p, namely, the set & of all the Hermitian operators in H that 
share with p the same set of Schmidt coefficients. Thus, we have that & p = & J % PI V(H). 

3. Characterization of the Schmidt equivalence classes 

Aim of the present section is to give a basic characterization of the Schmidt equivalence 
classes of states; see Definition [Q In this regard, one can adopt two different approaches. 
On one hand, one can try to characterize some 'natural' groups for which an action on the 
Schmidt equivalence classes is defined; see subsection 13. 1[ On the other hand, one can regard 
the Schmidt equivalence classes as manifolds and study their local geometry considering 
the action of the local orthogonal groups; see subsection 13.21 Although the results that we 
obtain are still somewhat 'preliminary', we think that it is worthwhile to report them since a 
description of Schmidt equivalence classes seems to be completely missing in the literature. 

3.1. Groups acting on the Schmidt equivalence classes 

Note that, since C U (H), £ R (7i A ) and £ R (7i B ) are real Hilbert spaces, the unitary 
(super)operators in these spaces belong to orthogonal groups. For instance, a unitary operator 
in C^(H) belongs to the orthogonal group 0(N A 2 N B 2 ). The class of unitary operators in 
£r(7~0 that are decomposable as the tensor product of two unitary operators in £ R (7i A ) and 
£r(^b), respectively — i.e. of the form f A ® % , with f A in 0(A r A 2 ), T B in 0(iV B 2 ) — will 
be denoted by O (g> O (H). It is clear that the maps in O <8> O (H) preserve the SC's of every 
element of (H), but, in general, T A <g> % (T>(TC)) <f. T>(TC). Note, moreover, that the set 
O <g> O (H) is a group (isomorphic to the direct product 0(iV A 2 ) x 0(iV B 2 )) with respect to 
the usual composition of maps. The orbit in C K {H), under the action of this group, passing 
through p E V(H) will be denoted by & p ; i.e. 

& p := {f(p) G £ n (H): T = T A ®f B 60®0 (H)}. (12) 
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Proposition 1 The orbit &p of the group O <8> 0{TL) through p e V{7i) coincides with the 
extended Schmidt equivalence class containing p: 

&p = &T- (!3) 

It follows that 

&p = &pnV(H). (14) 
Therefore, two states p and a in V(H) are Schmidt equivalent if and only if 

p = t A ®t B (&), (15) 
for some unitary operators T A and T B in C U (H A ) and C^(H B ), respectively. 

Proof: As already observed, given p e V(TL), for every T in the group O ® 0(H), T(p) 
belongs to Gf; hence, & p C Gf. On the other hand, let p = £)jj =1 h F£ <g> F fc B be a 
Schmidt decomposition of p and C = J2k=i G£ <g> Gl a Schmidt decomposition of an 
arbitrary element C of &J\ Then, for every couple of unitary operators T A and T B in C R (H A ) 
and C R (H B ), respectively^such that f A (F£) = G£ and T B (F fe B ) = G%, k = 1, . . . , d, we have: 
t A <g> t B (p) = (7. Hence, 6p D Sf. □ 

The set of all the maps of the form T A ®T B — with T A and T B unitary operators in Cj&(H A ) 
and C U (H B ), respectively — such that is stable under the action of T A <g> T B , i.e. such that 

t A ® t B (6^) C 6p, (16) 

is a semigroup (with respect to composition) with identity, contained in the group O ® O (H), 
semigroup which will be denoted by O (g> O st (p). The subset O (g> O inv (p) of O ® O st (p) 
defined by 

O ® O inv (p) := {t A ®f B GO® O st (p): (t A <g> t B )t GO® O st (p)} (17) 

is a group. It is easy to check that O (g> O inv (p) coincides with the subset of O ® O st (p) 
containing those maps that leave &p invariant, i.e. 

O <g> O inv (p) = {T A ®t B eO® O st (p): t A <g> t B (6^) = 6^}. (18) 

As an example of an operator belonging to O <g> Oi nv (p), for all p e T>(H), consider the linear 
map jf A ® J7" B : jCr(7Y) — > £ia(7i), with J" A and J B unitary operators defined by: 

J A . C R {H A ) 3A^ J A AJ A e £ R (H A ), J B . C R {H B ) 95h j B BJ B e C R (H B ), (19) 

where J A and J B are 'local' complex conjugations (i.e. selfadjoint antiunitary operators) in H A 
and H B , respectively. The maps J A and J B are partial transpositions, so that the map J A ® J B is 
the transposition associated with a tensor product basis in H A <S> H B (recall that transposition, 
as complex conjugation, is a basis-dependent notion). As it is well known, a transposition is a 
positive trace-preserving map (in short, PT map), and it is selfadjoint with respect to the HS 
scalar product. Therefore, the selfadjoint unitary operator J A ® J B is contained in the group 
O <g> O inv (p), for all peV(H). 

It is natural to wonder how states belonging to the same Schmidt equivalence class can 
be connected by physically realizable transformations. We will then consider the semigroup 
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with identity of PT maps in C R (TC), which will be denoted by PT(TC). It is worth defining 
the following subset of PT(7i:): 

PT e (7i) := {£ G PT(ft): £bijective, £~ l G PT(TC), 

8{&p) = & h Vp G V(H)}. (20) 

It is clear that the set PT q(TC) is a group. 

As already observed, unitary maps in C-^iTC) of the form T A <g> % — with %, and T B 
unitary operators in £ R (7i A ) and C R (TC B ), respectively — preserve the SC's of every element 
of Cn(TC), but, in general, T A ® T B (V(TC)) (jL V(H). It is therefore natural to consider the 
class of linear maps in Cjt(TC) belonging to the set O® O (TC) HPT(W). It is clear that this set 
is a semigroup (with respect to composition of maps) with identity contained in O ® O st (p). 
We will show now that it is actually a group which is a subgroup of O £g> O inv (p). We need 
two preliminary results; the proof of the first one is trivial. 

Lemma 1 The inverse of a bijective trace-preserving map from C R (TC) onto C R (TC) is trace- 
preserving. 

Lemma 2 A positive linear map T: C^(TC) — ► Cj&(TC) which is unitary transforms the convex 
cone of positive operators in H onto itself; therefore, is a positive map. Hence, in 
particular, if a linear map belonging to O ® O (TC) is a positive map, then its inverse is a 
positive map too. 

Proof: We will prove the statement by contradiction. Suppose that B = T(A) is positive, 
and assume that A is not a positive operator. Then, for some ip G TC, we should have that 
(ip,Aip) < 0; hence, ti(f(\^)(^\)f(A)) = tr(\ip){ip\ A) = (faAtJ;) < 0, where we 
have used the unitarity of T. On the other hand, since T is a positive injective linear map, 
^ = T(\ip){ip\) is a nonzero positive operator, so that it admits a decomposition of the form 
^ = Ylk=\ e k \4>k) {4>k\i where {4> k }^ =l is an orthonormal system and {e^}^ is a set of 
strictly positive numbers. Therefore, since B is positive, we have: 

K K 

tr^B) = Y,£ktr(\4 >k )( ( p k \B) = Y J £k(4>k,B4 >k ) >0. (21) 

k=l k=l 

But this is in contrast with the inequality tr(^ B) < previously found. □ 

From Lemma[T]and Lemma[2]one obtains immediately the following result: 

Proposition 2 The inverse of a linear map belonging to the setO®0 (Tt) fl PT (TC) belongs 
to this set too. Hence, O (g) O (TC) fl PT(TC) is a subgroup ofO® Oj nv (p). 

We now consider the set PTO ® PTO (TC) consisting of those linear maps in C R (TC) 
of the form £ A ® £ B , where £ A , £ B are linear maps in C R (TC A ), C R (TC B ), respectively, 
that are positive, trace-preserving and unitary. An analogous definition holds for the set 
CPTO <g> CPTO (TC), with the 'local maps' £ A , £ B assumed to be completely positive rather 
than simply positive. It is clear that the sets PTO ® PTO (TC) and CPTO <g> CPTO (TC) are 
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semigroups with identity. We will see that they are actually groups. Consider also the group 
of local unitary transformations 

0®0(H) := {f eO®0(H): f(A) = (U ®V)A(U ] ®V ] ),\/Ae C R (H), 
for some unitary ops. U, V in C R (H A ), C R (H B ) resp.}, (22) 

which is obviously a subgroup of both PTO® PTO (H) and OPTO® OPTO (H). In a similar 
way one defines the group of local unitary -antiunitary transformations O <E> O (TC) (include 
local antiunitary operators U, V in the r.h.s. of (1221) ). An example of a map that belongs to 
O <g> O (H), but not to O (g> 0(H), is the tensor product J A <g> J B of two partial transpositions 
J A and J B . 

Proposition 3 The set PTO ® PTO (TC) is a group (with respect to composition of maps); 
hence, it is a subgroup of the group O <E> O (H) fl PT(H). 

Proof: Given a map £ A ® £ B in PTO (g> PTO (H) — where £ A , £ B are linear maps in £ E (W fl ), 
£ R (7Y B ), respectively, that are positive, trace-preserving and unitary — the inverse map 
£ A l ® ^g 1 is positive; to show this, apply Lemma [Hand Lemma [2] identifying the (generic) 
Hilbert space H with the local Hilbert spaces H A and H B . □ 

A Kadison automorphism is a bijective map from V(K) — the convex set of density 
operators in a Hilbert space K — onto itself that is convex linear. We now state as a lemma a 
property of this kind of automorphisms that can be obtained 'by duality' from a well known 
result due to Kadison ll23l (concerning C* -algebras): 

Lemma 3 (Kadison) Every Kadison automorphism £: T>(K) — > T>(K,) is of the form 

£(p) = fpf\ VpeV(JC), (23) 
where T is a unitary or anti-unitary operator. 

Assume that the Hilbert space /C is finite-dimensional. Then, since any operator C G £r(X) 
can be written as C = C\ p\ — c 2 pi, for some p x , p 2 G T>(K) and some non-negative numbers 
ci, C2, it is clear that every Kadison automorphism £: T>(K) — > T)(K) extends (uniquely) in a 
natural way to a linear map in £ R (/C); conversely, a linear map in £ R (/C) which is bijective 
on V(K) can be regarded as a Kadison automorphism. 

Let us denote by O (TC) the group of unitary-antiunitary transformations in £j&(TC). We 
are now able to prove the following result: 

Theorem 2 The group 0<g>0(7^)nPT(7Y) coincides with the group 0®0(H)C) 0(H). The 
group PTO <8> PTO (H) coincides with the group O ® 0(H). The set OPTO ® OPTO (W) is 
a group which coincides with the group of local unitary transformations O <8> O (TC)- All the 
mentioned groups are subgroups ofPT & (H), and PT&(H) is a subgroup of 0(H). 

Proof: It is clear that the group O <g> O (H) 0(H) is a subgroup of O <g> O (H) n PT(W). 
On the other hand, by Lemma [Hand Lemma [2l a map in the group O <g> O (TC) fl PT(7Y) is a 
Kadison automorphism; hence, by Lemma[3l it is contained in 0®0 (H)dO (H). This proves 
the first assertion of the theorem. Next, given a map £ A <g) £ B in PTO g) PTO (7i) — where £ A , 
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£ B axe linear maps in Ar(7Y a ), £rCHb)> respectively, that are positive, trace-preserving and 
unitary — the maps S A and £ B are Kadison automorphisms. Hence, by Lemma[3l 

£ A (A) = UAU\ Vi e C U (H A ), 8 A (B) = VBV\ V5 g £ r (W b ), (24) 

for some unitary or antiunitary operators U, V in £r(7^ a ), C r (TC b ), respectively. Therefore, 
the group PTO®PTO (Tt) coincides with the group 0®0 (Tt). Note that, since an antiunitary 
operator is the composition of a unitary operator with a complex conjugation, if in (|24l) we 
let the operators U, V be antiunitary, we have that the maps £ A , £ B are the composition of 
unitary transformations with transpositions. Transpositions are positive but not completely 
positive maps; hence, the set CPTO ® CPTO (TC) coincides with the group of local unitary 
transformations 0®0(H). Finally, observe that the maps in the group PT 6 (H) are Kadison 
automorphisms. Thus, by Lemma[3l the last assertion of the theorem follows. □ 

3.2. Local analysis 

In Il24l . the authors considered the Schmidt decomposition of pure states, and, in that setting, 
analyzed the geometry of the sets of Schmidt equivalent pure states (that turn out to be 
differentiable manifolds). Since two pure states are Schmidt equivalent if and only if they 
are mutually convertible via local unitary transformations, they called such manifolds 'the 
manifolds of interconvertible states'. The aim of the present section is to apply the same 
line of reasoning to study the structure of the manifolds of Schmidt equivalent (mixed) states. 
Here we assume for simplicity that N A = N B = N, hence d = N 2 . This is not necessary for 
our purposes but it allows a convenient simplification of our formulae. 

As it is discussed below, there is a major limitation to the straightforward extension of the 
aforementioned results from pure to mixed states. The definition of SC's is based on the fact 
that density operators are elements of a vector space. On the other hand, they are constrained 
to be positive operators (of unit trace) because of their physical interpretation. This leads 
to the main difference between the results in [|24l and our forthcoming discussion: while 
in j[24| the geometry of the manifolds of interconvertible states can be characterized globally, 
here we are limited to a 'local' analysis (i.e. we are forced to consider transformations in a 
neighborhood of the identity). 

As stated in Proposition [H the transformations of the form p i— > T(p), where T belongs 
to O <E> O (TC), preserve the SC's. The converse is also true: if two operators have the same 
SC's then they are connected by a map in the group O (g> 0(H). Therefore, we can build 
the different equivalence classes acting locally on 'fiducial states' with the group O <8> O (TC). 
Actually, we need to consider the component connected to the identity of this group, which is 
isomorphic to (and will be identified with) the Lie group SO(d) x SO(d). We need then 
to impose two conditions on the local (i.e. close to the identity) transformations: a) the 
positivity of a fiducial state p must be preserved, and b) the trace of p has to be preserved 
as well. The first is a global constraint that does not 'reduce the number of dimensions', 
although characterizing a neighborhood of the identity in SO(d) x SO(d) preserving the 
positivity of a given fiducial state is a challenging open problem. The second constraint, 
normalization, is a local constraint that reduces the dimension of the manifold by one: 
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the condition that tr(p) = 1 amounts to fixing the projection of the vector p along the 
direction I. Let us consider an orthonormal basis in the real Hilbert space of the form 



{P <g> I B ,Ff <g> I B ,F 



PB PA 
1 3 ' 1 I 



*7}i<ij<d-i; clearly, tr(iP) = = tr(F B ), because 



of the orthogonality condition with I = I A <g> I B . Given a state p e V(H), we have: 

f A «P d_1 d_1 



Ft + E 



(25) 



i=l i=l i,j=l 

The infinitesimal action of the Lie group SO(d) x SO(d) on the basis elements is of the form 



JA,B 



is* + <t>T) f;> 



no 

JA,B 



+ e, A ' B F 



A, B 



(26) 
(27) 
1) 



where contraction of repeated indices is understood. Here A and B are real, (d — 1) x (d 
antisymmetric matrices and e A and e B are real (d — 1) -dimensional vectors. 

Hence, by applying an infinitesimal transformation in SO(d) x SO(d) with generators 
e B ) it is not difficult to see that the coefficient of I = fi A <8> I B — and hence tr(p) 



> A ,e A ; 



undergoes the following change: 

Str{p) = Ott + «fe B + ftrfej . 



(28) 



In this equation we have kept the leading non-trivial infinitesimal changes, discarding those of 
the type (e A ) 2 and (e B ) 2 . This is the correct expansion in the space of jets of regular functions 



of two vectors with nonvanishing gradients. Equation (1281) defines a surface similar to an 
hyperboloid. In order to visualize it one can simply take e A ' B = 5ue A ' B . For a simple choice 
of a, (3, the result is plotted in figure [Q Clearly, connection to the identity implies that only 
the branch containing the origin has to be considered. 




Figure 1. An example of the parameter space of allowed orthogonal transformations. Only 
the branch of the hyperboloid which contains the origin must be considered. 

Therefore, the vectors e A,B are constrained while the matrices </> A ' B are unconstrainedjjj 

§ We are neglecting here the case a A,B = 0, ftij ^ 0. In this atypical case we need to keep terms of O (e 2 ) 
in the transformation law of I A B . Once this is done one checks that the codimension is still one. The really 
degenerate case a = (3 = 0, instead, corresponds to the class of equivalence of a single point p = I A ® I B /d, 
with dimension 0. This point has to be considered as the 'tip', the extremal point of the space of states. 
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dp = — - 1 = d 2 - d - 1. (30) 



Assuming that the SC's are all different (i.e. that p is a 'typical state'), one finds out that the 
Schmidt equivalence class & p is locally diffeomorphic to 

SO(d) x SO(d) 

SO(2) ' K ] 

where the subgroup SO (2) is generated by the linear combination of the coefficients 
{ef , derived above. Then, the dimension of the manifold &^ is 

,d(d-l) 
2 

It is easy to see how this counting is consistent with the intuition that changing one of the 
SC's Ai, . . . , Ad of p brings out of the equivalence class &p. Indeed, subtracting the number 
of the SC's from the number of dimensions — (d 2 — 1) — d = d 2 — d — 1 — we obtain the 
dimensionality of a typical orbit. For example, take the case of two-qubits: N A = N B = 2. We 
have: d = 4 and d 2 — 1 = 15; hence, the manifold & p , for a typical bipartite state p E V(TC), 
is (15 — 4 = ll)-dimensional. 

If p £ V(H) is a non-typical state, the dimension of the manifold & p is smaller than 
d 2 — d — 1. In fact, in the case where the d SC's of p cluster into subsets of mi, ...,rrih 
identical values — by means of an argument analogous to the one adopted in ll24ll — one 
can check that there is a 'local stabilizer subgroup' isomorphic to SO (mi) x . . . x SO(m h ); 
therefore, in this case, & p is locally diffeomorphic to 

SO(d) x SO(d)/SO(2) 
SO(mi) x ... x SO(m fc )' 
Then, the dimension d p of & p is given, in general, by 

k=l 

In order to make the above argument rigorous, one should actually prove that the 
Schmidt equivalence classes are actually differentiable manifolds. In that case our previous 
argument would provide the correct dimension of such manifolds. As we learn from 
mathematicians [25 J, a standard tool for characterizing a subset of a differentiable manifold 
(like the manifold of Hermitian operators in TC) as a submanifold is the study of the (possible) 
Lie groups acting transitively on the given subset and of the associated stabilizer subgroups. 
Therefore, suitably improving the analysis of subsection 13.11 may allow to achieve such a 
remarkable result. 



4. Entanglement and symmetric polynomials in the Schmidt coefficients 

We will now consider the role played by Schmidt coefficients in the characterization of 
entanglement. Our starting point will be the realignment criterion (RC). It is natural to 
wonder if the whole set of the SC's of a bipartite state may allow a stronger characterization 
of entanglement with respect to the RC. 

As we have seen, the RC, like other separability criteria, is based on the evaluation of 
a single functional in terms of which a necessary condition for separability can be stated. 
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In the case of the RC, that functional coincides with the sum of the Schmidt coefficients. 
On the other hand the evaluation of a single functional might not be sufficient to determine 
completely the presence of entanglement j[22| (from a more general point of view, we may 
say, to characterize classical and quantum correlations). In particular, as the RC is only a 
necessary condition for separability, one is led to consider additional functionals in order to 
gain information about the presence of entanglement. 

Here we propose to consider the symmetric polynomials in the Schmidt coefficients, 
namely: 

m [i] = eLa 

m [/] = T, {h , i2 .., l} KK---K (33) 

M [d] = nLi x k- 

Notice that the RC involves the symmetric polynomial of degree one. 

A naive argument says that, if the sum of the Schmidt Coefficients is equal to S, their 
product is upper bounded by (S/d) d . Hence we have the following condition for a separable 
density matrix p: 

p separable M [d] < (^) , (34) 



which obviously defines a weaker separability criterion. Analogously, one can consider 
the symmetric polynomial of degree I and obtain the following necessary conditions for 
separability: 

p separable < y,(d) = (fj (~ \ . (35) 

In particular, if p has Schmidt rank R, we can write the conditions 

MU<( R )(±) 1 (36) 



J J V R , 

for / < R, while M [l] = for I > R. 

It is worth noticing that the symmetric polynomials {M'''} are in one-to-one 
correspondence with the Schmidt coefficients {Xk}- The inequalities (1351) are consequences 
of the RC. Hence, as separability criteria, they are weaker than the parent one. Section [6] will 
be devoted to the study of possible stronger generalization. 

Finally, we notice that in the approach followed in [12J, which makes use of the 
associated realigned matrix p R , the symmetric polynomials are the coefficients of the 
characteristic polynomial: 

d 

Xp n(x) = det(|p*| -xl) = J2 Mll] (\P R \)(-x) 6 ~ l (37) 

i=i 

where, with abuse of notation, we indicated with M^' (A) the principal minor of order I of the 
matrix A. 
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5. Quantum states and quantum maps 

This section is devoted to the application of the ideas presented in section @] to the study of 
quantum channels, i.e. completely positive trace-preserving (CPT) maps. In order to do that, 
we exploit the well known correspondence between quantum channels and quantum states 
(see ES 221 EH). 

Here we consider quantum systems with Hilbert spaces 7i A and 7i B , the set of states in 
the composite system V{H A £§) 7i B ), and the set of CPT maps from system B to system A, 
which is denoted CPT{H B , H A ). 

Given a CPT map £ G CPT(7i B , H A ), one can associate a state p G V(H A <g> H B ) in the 
following canonical way: 

£ — > p= (£®1)0), (38) 

where [3 = \<f>){<f>\ G T>(7i B ® H B ) denotes a maximally entangled state, for instance 
|0) = -J= J2a=i \ a ) an( ^ * s tne identical map in the system B. 

A CPT map £ is said to be entanglement breaking (EB) if £ ® X maps any state into a 
separable one E6l . One can show that a CPT is EB if and only if the map £ <g> I transforms a 
maximally entangled state into a separable one. It follows that the CPT map is EB if and only 
if the canonically associated state is separable. 

One can select a local orthogonal basis Fft^ = \m) (n\ and s = \(j){v\ and write the 
matrix elements of the CPT map in that basis as follows 

£(mn)(pu) = tr (Pfan)S ( F Qi V )\) ■ (39) 

It is easy to check that, with the canonical association (1381) . the matrix representation of 
the state p and the map £ are related in the following way: 

Following the definition in lfT2ll . it is immediate to recognize that, apart of the 
normalization factor 1/N B , the matrix expression of £ is identical to the realigned matrix p R 
(see equation ©). Identifying the CPT map with the realigned matrix of the corresponding 
density matrix, one can consider the characteristic polynomial 

P £ {x) = det(|£|-xl), (41) 

where \£\ = \f£X£ and £ = £{ij)( a p)- To fix the ideas, let us consider the case in which £ has 
full rank. One obtains from (|34|) the following necessary condition for £ to be entanglement 
breaking: 

det(|£|)< f^j . (42) 

The determinant of quantum channels was also considered in [27 J, in which some of 
its properties were presented and discussed in the context of factorization of CPT maps. The 
present result relates a geometric property of the map, such as the rate of contraction of volume 
(which is equal to det (|£|)) to the property of being entanglement breaking. Analogously, 
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from (1361 ), if the matrix £(jjv a( g) has rank R, we can write the following necessary conditions 
for E to be entanglement breaking: 

M [l H\£\)<( R )(^) 1 (43) 

for I = 1 . . . R. 



I R 



6. Beyond the realignment criterion 

In section |4] we introduced a family of separability conditions which are weaker than (or 
equivalent to) the RC. In this section we argue about the possibility of extending the family of 
separability criteria defined in ((351) in order to write criteria which are independent of the RC. 
As a first step in this direction, we may ask whether it is possible to find strict upper bounds 
xi(d, D) such that: 

p separable M«<^(d,D)< Q . (44) 

However, in the following we consider a weaker statemenjjj, namely 

MW < 1 < f,(d, D) < (fj (^) (45) 

which establishes a strict upper bound for the functionals over the set of states satisfying 
the RC. 

The following proposition holds true: 

Proposition 4 The upper bounds xi(d, D) in equation rf?5l) exist for D < d 3 . 

Proof: The proposition is proven by contradiction. Let us suppose the existence of a density 
matrix p such that M' 1 ' < 1 and the inequalities in (1351) are saturated. That implies that 
the density matrix has maximum rank, R = d, and all its SC's are all equal to d _1 . Hence, 
referring to the singular value decomposition in equation (flOl) . we can write the Schmidt 
decomposition of the density matrix in the following way: 

1 d 

P0(ai)(l3j) = "J U (af3)(ld) V \kl)(ij)' ( 46 ) 

k,l=l 

where v (ij)(i'j') and u^p)^ p) are respectively the entries of the unitary matrices V and U (see 
equation (flOl)), with dimension iV A 2 and iV B 2 . We introduce the following notation 

U( a /3)(M)V* (kl){ij) = (Vij^Uap), (47) 

where u a p and indicate the vectors respectively defined as the rows and the columns of the 
matrices U and V. We have: 

tr (Po) = -^^2(v»,u aa ) = -^(^2v ih ^2u aa ), (48) 

i,a i=l a=l 

|| Notice that d45l ) is weaker as separability criterion, while it is stronger in the sense that d45l l implies (l44l) . In 
particular x;(d, D) < i;(d, D) 
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Figure 2. A pictorial representation of the set of the equivalence classes, for the generic case 
d 7^ D (a); and for the degenerate one d = D (b). The Schmidt coefficients are represented 
along the axes. States are contained in the colored region, which fulfills the constraints on the 
purity ©. The region on the bottom-left of the solid line fulfill the RC (hence separable states 
are contained in this region). The region on the bottom-left of the dashed line contains the 
states that fulfill one of the naive inequalities d35l l for some Z > 1. Finally, the region on the 
bottom-left of the dot-dashed line contains the states that fulfill the proposed criterion RCL 
Notice that the RC (or RC1) bounds separable states with a hyperplane, while the higher order 
RCZ, for I > 1, bounds the set of separable states with non-linear hypersurfaces. Notice that 
the RCTs can in principle allow a finer characterization of entangled states. 



and we obtain: 

tr {po ) = | tr ( fi0 ) | = KSijl^MI < [SiMfij^ = 3^, (49 ) 

(the last equality holds true since u aa and va are two systems of orthonormal vectors) which, 
for D < d 3 , is in contradiction with the hypothesis that p has unit trace. Since the set of states 
satisfying the RC is compact and the symmetric polynomials are continuous functionals, 
the lower upper bounds x;(d, D) do exist. □ 

In the following we indicate with RC/ the suggested criterion 

p separable =>■ M V] < a?,(d, D). (50) 

Figure [2] shows a pictorial representation of the relation between the parent RC, the weaker 
criteria (1351) and the proposed extensions RC/. As the figure suggests, the RC/ can in principle 
be used, with respect to the information given by the RC, as refinements of the knowledge 
about the region of separable states. 

It is worth noticing, however, that there are still two main open problems: 

(i) the actual values of the upper bounds x/(d, D) (as well as xi(d, D)) are still undetermined; 

(ii) it is not clear whether the criteria RC/ are independent of the RC, i.e. if there are 
entangled states such that RC is not violated while RC/ is for some / > 1. That is 
equivalent to the strict inequality x;(d, D) < Xi(d, D). 

In the next section we face these problems with a numerical approach. We are going to restrict 
our discussion to the case of lower dimensional systems, namely 7i A ®7f B = C 2 ®C 2 , <C 2 ®C 3 . 
For this cases one can exploit the fact that the PPT (positive partial transpose) criterion j8]| is 
necessary and sufficient for separability. 
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7. Examples for low dimensional systems 

In analogy to what can be done for the RC (see for instance lfT9lO . one could determine the 
value of the strict upper bounds Xi(d, D) by convex linearity starting from the properties of 
pure separable states. Nevertheless, it is worth noticing that that can be a rather difficult task 
since the symmetric polynomials M^' are not easy to manipulate with respect to the convex 
structure of the set of separable states. For this reason, in the following we present a numerical 
analysis which allows to present some interesting results. 

For a preliminary analysis of the potentialities of the proposed family of criteria, they 
have been numerically tested in the case of a bipartite qubit-qubit and qubit-qutrit system. 
A numerical search of the upper bounds can be done exploiting the PPT criterion [8J. The 
constraints p > and p TA > 0, where Ta indicates the partial transposition, are known to be 
necessary and sufficient to characterize separable states in the low dimensional cases. Hence, 
the determination of the lower upper bounds xi(d, D) reduces to a problem of constrained 
maximization. 

We have numerically estimated the maxima of the functions over separable states 
(hence determining estimates for xi(d, D)), and over the set of states satisfying the RC (hence 
estimating the upper bounds xi(d, D)). For a qubit-qubit system (d = D = 4) the results are 
shown in table Q] together with naive bounds in (1351) . The analogous quantities are shown in 
table [2] for the case of a qubit-qutrit system (d = 4, D = 9). In the latter case we found 
xi(d, D) < xi(d, D) suggesting that the criteria RC/ can be in principle stronger than the RC. 





I = 1 


I = 2 


/ = 3 


1 = 4 




1 


0.3750 


0.06250 


0.003906 


x,(4,4) 


1 


0.3333 


0.04630 


0.00231 


x,(4,4) 


1 


0.3333 


0.04630 


0.00231 



Table 1. For a qubit-qubit system, the table shows the upper bounds on the sums of the 
symmetric polynomials M'''. yz(4) denotes the naive bounds d35l l. x&(4, 4) denotes the 
numerically estimated strict bounds over the states satisfying the RC d45l l. Finally, x;(4, 4) 
indicates the numerically estimated strict bounds over the set of separable states (l44l . Notice 
that in this case we found x;(d, D) = xi(d, D) suggesting that the criteria RCl cannot be 
independent of the RC. 

Figures [3] and |4] show the maximum of the functionals for I = 2, 3, 4, computed for fixed 
values of M^, respectively for qubit-qubit and qubit-qutrit system, as functions of the value 
of M [1] . The maxima are computed over generic states and over separable states. In the case 
of qubit-qutrit system, the plots show the region in which the criteria RC/ can in principle be 
stronger than — or independent to — the RC. 

To conclude this section, we consider the case of two-qubit (generalized) Werner states, 
of the form: 



P P = p|0X0l + (l-#/4, 



(51) 
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/ = 1 


I = 2 


/ = 3 


1 = 4 


yi(4) 


1 


0.3750 


0.06250 


0.003906 




1 


0.3583 


0.05533 


0.003133 


^(4,9) 


1 


0.3469 


0.05249 


0.00291 



Table 2. For a qubit-qutrit system, the table shows the upper bounds on the sums of the 
symmetric polynomials M^'. y;(4) denotes the naive bounds d35l >. 5j(4, 9) denotes the 
numerically estimated strict bounds over the states satisfying the RC d45l l. Finally, xi(4, 9) 
indicates the numerically estimated strict bounds over the set of separable states (|44| |. Notice 
that in this case we found xi(d, D) < xi(d, D) suggesting that the criteria RC/ can be in 
principle independent of the RC. 




Figure 3. Numerically estimated upper bounds of the symmetric polynomials (a), 
(b), and M' 4 ' (c), as functions of the value of m' 1 ' for a qubit-qubit system: naive upper bounds 
(dash-dotted line), numerically estimated upper bounds over the set of separable states (solid 
line) which coincides with upper bounds over the set of all states (separable and entangled). 




Figure 4. Numerically estimated upper bounds of the symmetric polynomials (a), 
(b), andM' 4 ' (c), as functions of the value of for a qubit-qutrit system: naive upper bounds 
(dash-dotted line), numerically estimated upper bounds over the set of separable states (solid 
line) and over the set of all states (dashed line). Notice the presence of a region in which the 
functionals have a lower upper bound over the set of separable states. 
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for p G [0, 1], where \ (f>) indicates a maximally entangled pure state. The Schmidt coefficients 
of the state (I5TI) are easily calculated tobe{l/2,p/2,p/2,p/2}, yielding 

M [1] = (l+3p)/2 
M [2] = 3(p + p 2 )/A 

= (3j9 2 +j9 3 )/8 K } 

M [4] = p 3 /16. 

Notice that the symmetric polynomials are monotonically increasing functions of the state 
parameter p. It was shown in [|28l that the realignment criterion is necessary and sufficient for 
this family of states (indeed it is so for all the two-qubit states with maximally disordered 
subsystems). The state in (I5TI) is known to be separable for p E [0, 1/3] and entangled 
otherwise. We can compute the maximal value of the symmetric polynomials over the 
separable states in that family. These maxima are reached in correspondence of the value 
p = 1/3, hence yielding: 

M [1] < 1 
M [2] < 1 /3 

p p separable =► ^ ^ ^ ^ Q QmQ (53) 

M [4] < 1/432 ~ 0.00231. 

Notice that the values in (l53l) computed for two-qubit generalized Werner states saturates the 
numerical estimated upper bounds reported in the tabled] 



8. Conclusions 



The main goal of the present paper is to bring attention to the Schmidt coefficients of a 
bipartite density operator, and to their role for entanglement detection. The notion of Schmidt 
equivalence classes has been introduced and a preliminary characterization of such classes 
has been provided. 

We have presented a family of separability criteria, written in terms of the Schmidt 
coefficients, which are derived from the realignment criterion. These separability criteria 
are consequence of the fact that the symmetric polynomials in the Schmidt coefficients are 
upper bounded on the set of separable states. 

The application of that family of criteria to the study of quantum channels determines a 
relation between a physical feature, such as the preservation of entanglement under the action 
of the channel, and a geometrical quantity, such as the determinant — or the sum of principal 
minors of order I — of a corresponding matrix. 

We conjecture — also with support of numerical examples — that a strengthened version 
of these criteria, independent of the realignment criterion, exists. In particular, we have given 
numerical examples for the case of the qubit-qutrit system. These numerical results are of 
course not sufficient for achieving independent separability criteria. However, they can open 
the way to an analytical determination of stricter upper bounds on the symmetric polynomials, 
and this may eventually lead to new separability criteria. 
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